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GROUND  WAVE  PROPAGATION  IN  THE  PRESENCE  OF  SMOOTH  HILLS  AND  DEPRESSIONS* 


L.  B.  FelsenandA.  Green 
POLYTECHNIC  INSTITUTE  OF  NEW  YORK 
Farmingdale,  N.  Y.  11735  USA 


I.  INTRODUCTION 


Electromagnetic  waves  propagating  along  the  ground-air  interface  are  influenced  markedly  by 
terrain  features.  Previous  studies  have  dealt  extensively  with  abrupt  changes  in  surface  properties 
(modeled  in  terms  of  surface  impedance),  or  abrupt  changes  in  height  as  represented  by  precipitous 
hills  or  valleys,  on  the  convexly  curved  earth's  surface.  When  the  terrain  changes  gradually,  as  in  the 
presence  of  smooth  hills  or  depressions,  (Fig.  1),  it  is  necessary  to  consider  the  effects  of  variable 
concave  surface  curvature,  and  more  generally  those  of  a gradual  changeover  from  concave  to  convex, 
on  the  signal  propagation  characteristics.  Of  special  interest  are  the  fields  observed  on  or  near  the 
surface  when  the  source  is  also  situated  on  or  near  the  surface.  This  paper  is  concerned  with  a review 
of  earlier  work  on  the  concave  surface  problem,  and  with  a summary  of  our  current  effort  to  adapt  and 
extend  these  results  to  ground  wave  propagation. 

The  most  fundamental  difference  between  the  excitation  of  concave  and  convex  surfaces  of  large 
(compared  to  the  wavelength)  radius  of  curvature  is  the  absence  in  the  former  of  a geometrical  shadow 
region,  from  which  the  source  is  invisible.  This  circumstance  gives  rise  not  only  to  a more  intricate 
geometric -optical  field  comprising  multiply  reflected  rays  (Fig.  2)  but  also,  in  an  alternative  guided 
wave  description,  to  the  presence  of  whispering  gallery  modes  which  cling  to  the  surface  (Fig.  3)  and, 
in  the  absence  of  dissipation,  experience  no  attenuation.  Their  counterparts  on  a convex  surface,  the 
creeping  waves,  lose  energy  by  tangential  shedding  along  the  propagation  path.  The  problem  is  compli- 
cated further  by  the  fact  that  ray  optics  is  incapable  of  providing  the  field  solution  since  the  caustics  of 
multiply  reflected  rays  pile  up  near  the  boundary  (Fig.  2)  and  thus  invalidate  the  geometric-optical  field 
evaluation  there.  A field  representation  in  terms  of  whispering  gallery  modes  only  (with  inclusion  of  a 
continuous  mode  spectrum  for  some  surfaces),  while  valid,  is  inconvenient  for  calculation  for  large 
separation  of  source  and  observation  points  on  a large-radius  surface  since  the  number  of  modes  required 
can  be  substantial.  These  problems  do  not  arise  on  a convex  surface  where  the  distant  field  in  the  shadow 
region  is  represented  compactly  by  the  dominant  creeping  wave. 

To  gain  a better  physical  as  well  as  quantitative  understanding  of  these  aspects  of  wave  propagation 
on  a concave  surface,  intensive  studies  have  been  carried  out  on  the  simplest  prototype  configuration,  the 
interior  of  a perfectly  reflecting  circular  cylinder  (KINBER,  B.  YE,  1961;  BABICH,  V.  M.  and  BULDYREV, 
V.  S.  , 1 972;  WASYLKIWSKYJ,  W.  , 1975).  These  two-dimensional  field  calculations,  corresponding  to 
excitation  by  an  axial  line  source,  are  reviewed  in  Section  II.  A peculiarity  of  the  cylindrical  geometry 
is  the  presence,  in  addition  to  the  whispering  gallery  modes,  of  a continuous  guided  mode  spectrum 
which  arises  because  of  spurious  reflections  from  the  radial  coordinate  origin.  Elimination  of  these 
spurious  contributions  leads  to  an  asymptotic  field  representation  in  terms  of  an  integral  which  can  be 
manipulated  so  as  to  exhibit  ray-optical  contributions,  whispering  gallery  mode  contributions,  a mixture 
of  these,  or  a formulation  containing  a reduced  canonical  integral  analogous  to  the  Fock  integral  for  con- 
vex surfaces.  The  most  effective  choice  depends  on  the  parameters  of  the  problem.  The  results  so 
obtained  can  be  generalized  to  apply  to  arbitrary  concave  surface  shapes  provided  that  the  radius  of  cur- 
vature changes  slowly  over  a wavelength  interval.  The  validity  of  these  variable -curvature  solutions  has 
been  verified  by  comparison  with  exact  calculations  performed  for  a parabolic  contour  (BULDYREV,  V.  M. 
and  LANIN,  A.  I. , 1975).  These  earlier  investigations  are  reviewed,  relevant  numerical  results  pre- 
sented and  their  physical  implications  discussed. 

We  have  extended  the  analysis  to  track  the  field  behavior  when  the  radius  of  curvature  of  a circular 
cylindrical  surface  between  fixed  source  and  observation  points  changes  from  concave  to  convex;  this 
transition  from  a curved  to  a flat  boundary  and  beyond  was  not  included  in  the  earlier  studies.  We  have 
also  introudced  the  generalization  for  a lossy  ground  modeled  by  a surface  impedance  boundary  condition 
(Section  UI). 

The  two-dimensional  results  are  useful  for  the  investigation  of  fields  with  broad  phase  fronts  that 
do  not  vary  appreciably  along  the  axial  direction  of  an  elongated  depression  (gulley)  or  smooth  ridge. 

When  the  incident  phase  front  is  strongly  curved,  as  is  the  case  when  a dipole  radiator  is  located  within 
or  near  the  ground  perturbation,  it  is  necessary  to  calculate  the  three-dimensional  fields.  These  modifi- 
cations are  under  study. 

Still  to  be  explored  is  the  field  description  along  surface  contours  that  vary  smoothly  from  concave 
to  convex.  This  important  problem  is  presently  being  investigated,  as  is  the  extension  of  the  analysis 
to  three-dimensional  ground  perturbations  exemplified  by  spherical  depressions. 


* 


This  work  was  performed  under  Contract  No.  DAHC  04-75-G=01  52,  U.  S.  Army  Research  Office  - 
Durham,  N.  C. 
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U,  PERFECTLY  CONDUCTING  CIRCULAR  BOUNDARY  (TWO-DIMENSIONAL) 

A.  Green's  Function  Formulation 

We  summarize  here  the  results  obtained  in  the  literature  on  the  prototype  problem  of  line  source 
radiation  inside  a perfectly  conducting  circular  cylinder  of  radius  a.  To  make  the  analysis  relevant  to 
the  propagation  along  concave  surfaces  in  general*  it  is  necessary  to  remove  the  azimuthal  periodicity 
imposed  on  the  fields  in  the  cylindrical  g = (p,  0)  geometry,  where  p represents  the  radial  and  0 the 
angular  (azimuthal)  coordinate.  This  can  be  done  by  extending  the  range  of  the  6 coordinate  from  its 
physical  periodic  domain  0 < 0 < 2Tr  into  an  infinite  domain  -oc  < 0 < oo  (FELSEN,  L,  B.  and  MARCUVITZ, 
N,  , 197  3).  Such  an  extension  implies  that  angularly  propagating  waves  originating  at  the  source  point 
(p',0*)  are  outgoing  toward  101  = oo;  this  "angular  radiation  condition"  can  be  realized  by  placing  at 
some  angular  location  away  from  the  source  angle  0*  a perfect  absorber  for  angularly  propagating 
waves.  Such  an  absorber  has,  however,  the  undesired  property  of  generating  diffraction  at  the  radial 
coordinate  origin  p = 0.  Therefore,  when  considering  propagation  phenomena  ascribable  only  to  the 
cylinder  surface,  it  is  desirable  to  remove  the  spurious  diffraction  effects  from  the  total  field  solution. 

The  line  source  Green's  function  G(jg,  £')  in  the  cylindrical  domain  0 < p < a,  -co  < 0 < oo,  can  be 
constructed  by  the  method  of  separation  of  variables  and  expressed  in  its  most  general  form  as  a contour 
integral  involving  the  two  one-dimensional  characteristic  Green's  functions  gp  and  g^  for  the  radial  and 
angular  domains,  respectively  (FELSEN,  L.  B.  and  MARCUVITZ,  N, , 1 973),  As  noted  in  Section  I,  when 
the  Green's  function  is  represented  in  terms  of  angularly  propagating  waves,  the  corresponding  eigenmode 
spectrum  in  the  radial  domain  involves  a discrete  (whispering  gallery  mode)  and  continuous  portion.  The 
latter  accounts  in  part  for  propagation  phenomena  associated  with  the  portion  of  the  concave  boundary 
lying  between  the  source  and  observation  points  (this  is  desired  contribution),  and  in  part  for  the  spurious 
diffraction  effects  of  the  angular  absorber.  When  the  spurious  effects  are  removed  and  cognizance  is 
taken  of  the  high-frequency  nature  of  the  analysis,  one  may  show  that  the  relevant  propagation  character- 
istics along  the  portion  of  the  concave  boundary  lying  between  the  source  point  Q and  the  observation 
point  P in  Fig.  2 are  contained  in  the  partial  Green's  function  (BABICH,  V,  M,  and  BULDYREV,  F.  S. , 
1972) 


G„(e.  £') 


1 I 

i(iTka)^  *^*(ka)  P (ka) 


(1) 


where  k is  the  free-space  wavenumber,  the  prime  on  the  cylinder  functions  denotes  the  derivative  with 
respect  to  the  argument,  and  a time  factor  exp(-iu>t)  is  implied.  It  has  been  assumed  in  (1)  that  P and 
Q both  lie  on  the  boundary  (i.  e. , p - p'  = a)  and  that  the  boundary  condition  requires  that  the  normal 
derivative  of  the  Green's  function  vanishes  at  p = a;  this  makes  Gq  proportional  to  the  axial  component 
of  magnetic  field.  The  contour  C and  the  singularities  of  the  integrand  in  the  complex  i/  -plane  are 
shown  in  Figure  4.  Contributions  from  the  pole  singularities  arising  from 


J'  (ka)  = 0 , m = 1,  2,  . . . (2) 


are  found  to  describe  whispering  gallery  modes.  While  (2)  has  an  infinite  number  of  real  solutions  as 
indicated  in  Fig.  4,  only  those  with  Re  > 0 represent  spectral  contributions  in  the  angular  trans- 
mission representation  which  includes  also  a continuous  opectrum. 

B.  Rav-optical  Representation 

The  integralin  (1)  may  be  dealt  with  in  various  ways  which  exhibit  the  role  of  whispering  gallery 
modes,  ray-optical  fields,  or  a combination  of  these,  as  noted  in  Section  I.  If  an  evaluation  is  sought 
by  the  method  of  steepest  descent,  one  replaces  and  = (H^  (2))/  2 by  their  Debye 

asymptotic  forms 


J ^ Jcos  w - - w)  sinwj  + i 

Re  w > 0,  I ka  fcos  w - (-j  - w)  sin  wj  I ^ > 1 


(3) 

(3a) 


where  t/  = ka  sinw.  Utilizing  the  traveling  wave  expansion 

H.  (1) 


1 


2 


r (ka) 


H 


. (2) 


(ka) 


.1 


(ka) 


(-1)*' 


H’  ^^’(ka) 

v 


|r|<l  , 


(4) 


in  conjunction  with  (3),  one  reduces  (1 ) to  a series  of  integrals 
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G -- 
o 


00 

n>0  C 


where 


q^(w)  = |<()-0'|  sinw  + 2(n+l ) Jcos  w - (^  - w)  sin wj 


The  saddle  points  of  q (w),  as  obtained  from  dq  / dw  = 0,  are 


V = " . k -0*1 
sn  2 2(n+l) 


— ir 
*8  ■ 2 


(5) 


(6) 


(7) 


and  thus  lie  on  the  real  w axis  between  w = 0 and  w = ir/2.  A typical  path  through  the  nth  saddle 
point  is  shown  in  Figure  5.  The  original  path  mapped  from  the  i/  plane  can  be  deformed  into  C . The 
Debye  approximations  in  (3)  are  valid  in  this  relevant  region  of  the  complex  w-plane.  Use  of  the  conven- 
tional saddle  point  formula  for  evaluation  of  the  integral  in  (5)  then  yields  the  geometrical  ray  series 


G 

o 


ikD 

n 


where 


D^  = 2(n+t ) a sin  ( k - 0'  | / 2(n+l )) 


(8) 


(8a) 


each  term  of  which  corresponds  to  a direct  or  multiply  reflected  ray  as  depicted  in  Figure  2(a).  The 
result  in  (8)  could  have  been  constructed  directly  by  ray-optical  techniques. 

C.  Mixed  Rav-optical  and  Whispering  Gallery  Mode  Representation 

Although  the  series  in  (8)  formally  contains  ray  contributions  with  an  arbitrarily  large  number  of 
reflections,  these  are  suspect  since  the  saddle  points,  from  which  they  are  derived  all  cluster  about 
w = it/  2.  Thus,  the  asymptotic  method,  whereby  each  saddle  point  is  treated  as  isolated,  is  inadequate. 
Moreover,  (3)  becomes  invalid  as  w-»tr/ 2,  thereby  invalidating  the  simplification  of  the  integrand,  on 
which  the  saddle  point  evaluation  is  based.  It  is  therefore  necessary  to  truncate  the  number  of  legitimate 
ray-optical  terms  at  some  n = N such  that  Wg^  is  sufficiently  less  than  ir/  2.  With  this  in  mind,  one 
employs  instead  of  (4)  the  partial  expansion 


1 

J'  (ka)  ■ 
1/ 


H' 

1/ 


N 


Z 


n=0 


J'  (ka) 
1/ 


(-1) 


N+1  N+1 


(9) 


When  (9)  is  substituted  into  (1 ), 


one  may  write 


G = G 1^  + 
o oN  N 


(10) 


where  is  given  by  (5)  with  the  series  truncated  at  n = N,  and 


. ^ikak-0'|si..w 

: J (21 *■  (-1)  coswdw 

C*  H'  ''‘'(ka)J' 


.(ka) 


(ID 


where  v = ka  sinw.  The  asymptotic  evaluation  of  Cpw  evidently  yields  ray  contributions  having  expe- 
rienced up  to  N reflections.  The  remainder  integral  incorporates  the  cumulative  effects  of  ray 
fields  having  been  reflected  more  th^  N times.  It  should  be  noted  that  when  the  observation  point 
approaches  the  source  point  so  that  10  -0'|-»  0,  one  has  Wg-»B/2  even  for  n * 0 and  the  ray  repre- 
sentation fails  altogether. 

Several  options  are  open  for  dealing  with  Rj^.  One  possibility  is  to  deform  the  integration  path  so 
as  to  minimise  |R{,(|.  If  C'  is  deformed  into  in  Fig.  5,  M residues  at  the  poles  Vm  of  the  integrand 
lying  between  C'  and  Ciq  must  be  extracted.  Thus,  noting  that  r = -1  at  (see  (2)  and  (4)), 
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G . . + 
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MN 


(12) 


where 


°oM  * ka 


,u  -«■ 


m»l 


(13) 


with  v—s  ka  aln  w,^  and  being  the  same  as  except  for  the  replacement  ofC  by  C^.  Each  term  in 

the  series  represents  pre^sely  one  of  the  whispering  gallery  modes  which  make  up  the  discrete  spectrum 
in  the  previously  mentioned  Green's  function  representation  in  terms  of  angularly  propagating  waves. 
Since  w^<  ir/2.  one  may  employ  the  Debye  approximation  in  (3)  to  simplify  (12).  Thus,  from  (10)  and  (12), 


G 

o 


G w + 
oN 


G + 

oM 


MN 


(14) 


D.  Whispering  Gallery  Mode  and  Continuous  Spectrum  Representation 

The  spectral  representation  comprising  the  total  discrete  (wliiapering  gallery  mode)  portion  and  a 
continuous  portion  may  be  obtained  directly  from  (t ). 

As  remarked  in  Section  A,  only  those  pole  singularities  liaving  > 0 yield  spectral  contributions 
in  the  angular  transmission  representation.  When  C in  Fig.  4 is  deformed  into  a contour  extending 
along  the  imaginary  v axis,  all  necessary  pole  singularities  are  picked  up  and  one  is  left  with 


G 
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G w + 
oM 


_L 


i(n  ka) 


2 


tigr 


,i‘^l»  •»'ldv 

H*  ^^Nka)  J*  (ka) 


(15a) 


where  Gq]^  is  given  by  (1  3)  and  M is  the  total  number  of  poles  satisfying  Re  > 0,  The  second 
member  of  (1  5s)  is  the  continuous  spectrum  mentioned  in  Section  A except  that  the  spurious  diffraction 
effects  of  the  angular  absorber  are  not  present.  By  methods  similar  to  Wasylkiwskyj  (WASYLKIWSKYJ, 
W.  , 1975),  this  integral  may  be  simplified  considerably  to 


i J dv 


(I5b) 


which  is  valid  for  |0  -0'|  arbitrarily  small.  It  is  also  pointed  out  that  (l5b)  may  be  written  as 


i [H^(kaU-0'I)  - Y^(ka|,/)-0'I  ) ] (l5c) 


where  H_  is  the  Struve  function  and  Y the  Neumann  function  which  provides  the  correct  singular 
behavior  as  0 -»  0'  . 

E.  Optimum  Mix  of  Rays  and  Modes 

The  utility  of  the  mixed  ray  - mode  representation  in  (14)  depends  on  the  best  choice  of  N and 
hence  M so  as  to  minimize  Rvjw  and  render  it  negligible.  Utilizing  the  asymptotic  approximations  in 
(3)  to  simplify  the  integrand  of  Rj^,  Wasylkiwskyj  (WASYLKIWSKYJ,  W.  , 1975)  has  concluded  that 
^MN  neglected  at  observation  points  characterized  by  Wgj^j  = (wj^^+wj^+j )/ 2,  i.  e. , for  saddle 

points  well  removed  from  the  pole  singularities.  The  justification  is  fairly  loose  and  requires  further 
study.  It  is  likely  that  a precise  choice  is  unimportant  for  very  large  ka  and  large  separation  of  source 
and  observation  points.  However,  for  moderate  values  of  ka,  or  when  the  observer  approaches  the 
source,  the  problem  must  be  addressed  in  a different  manner. 

Returning  to  the  original  expression  for  Gq  in  (1),  one  may  explore  under  what  conditions  the 
contributing  range  of  the  integrand  is  localized  so  that  approximation  methods  may  be  used  effectively. 
To  this  end,  the  cylinder  functions  are  first  replaced  by  their  uniform  asymptotic  approximations  in 
terms  of  Airy  functions; 


H'  *^’(ka)  - iBi'(-TV/^)1 


n=  ka 


(16) 


J* 

V 


(ka)  - 


2 cos  w 


T 


(17) 
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- (cos  w - - w ) sinw) 


1/3 


Re  w > 0 


2 2/3 

When  T (ka)  ' > >1,  these  representations  reduce  to  those  in  (3).  Examining  the  behavior  of  the 

integrand  over  the  contour  C in  the  i/-plane  or  the  contour  C in  the  w-plane,  Babich  and  Buldyrev  have 
shown  (BABICH,  V.  M.  and  BULDYREV,  V.  S. , 1972)  that  the  principal  contribution  to  the  integral  arises 
from  the  vicinity  of  v = ka  or  w = 'ir/2  provided  that  the  arc  length  parameter 


s = a|0  - 0’  I 


(18) 


is  bounded  by  the  inequality 


♦ 


_2 


<Y<2./r  (^) 


6 


(19) 


where  ^ is  a quantity  of  0(t)  for  moderate  ka  and  6 is  small  and  positive.  Then  changing  variables  to 
t via 


and  using 

H'<2’(ka)--  (j^)  w-(t) 


o 2/3 

J^**^**--  Ai'(t) 

where 


(20) 


(20a) 


W2(t)  = Ai(t)  + i Bi(t) 
one  may  write  in  the  form 


G 

o 


iks 
e 

2iir  ^ka 


(— ) 
' 2 ' 


2/3 


J 


w^(t)Ai'(t) 


dt 


(20b) 


(21) 


where  the  contour  C in  the  t-plane  is  inferred  from  (20)  and  Fig,  4 by  observing  that  v * ka  corresponds 
to  t B 0 (note  that  the  mappings  and  tp  of  and  Vn,  respectively,  signify  the  zeros  of  Ai’(t)  and 
w^(t).  Extracting  ray  optical  contributions  as  in  (10),  when  relevant,  the  form  for  Rjj  analogous  to  (11) 
becomes 


«N  = 


iks 

e 

2iit  ^ka 


<¥> 


2/3 


J 

C 


Jjll 


N+1 


Ai'(t)  w^t) 


h<^>i 

[wyt)J 


N+l 


iVt 


dt 


(22) 


Buldyrev  and  Lanin  state  as  a criterion  for  the  inclusion  or  not  of  the  ray-optical  series  G.u  the 
inequality 


-1  < (N+l)  < 


(23) 


with  6 taken  as  (l/21).  When  the  left-hand  side  of  the  inequality  is  negative,  no  ray-optical  contributions 
are  included.  The  canonical  integral  Rja  in  (22)  has  been  tabulated  for  N > -1,  0,  2 (BULDYREV,  V.  M. 
and  LANIN,  A.  I.  , 1975). 

* Ref.  2 actually  deals  with  the  Dirichlet  problem  ■ 0 on  the  boundary.  The  validity  of  exactly  the 

same  condition  for  the  Neumann  function  used  here  remains  to  be  checked. 
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F.  Near  Field  and  Infinite  Plane  Limita 

The  preceding  results  cannot  account  for  the  limiting  case  y-^  0 in  (18),  which  arises  either  when 
the  observation  point  approaches  the  source  point  (s 0 with  a fixed)  or  when  the  radius  of  curvature 
becomes  arbitrarily  large  (a  « with  s arbitrary  but  fixed).  The  latter  limit  traces  the  transition 
from  a concave  to  a plane  surface  and,  if  a is  allowed  to  become  negative,  from  concave  to  convex.  We 
have  derived  the  transition  functions  for  these  cases. 

We  begin  with  Gq  in  (21)  and  use  the  Wronski  ,n  relation  for  the  Airy  functions  together  with 
Cauchy's  theorem  applied  to  the  upper  half  of  the  v-plane  to  infer  that,  equivlently. 


with  the  contour  C passing  below  the  real  t axis.  Regarding  the  integral  in  (24)  as  a Laplacg  integral, 
the  principal  contribution  for  small  y arises  from  the  range  of  large  t.  Shifting  the  contour  C so  that 
this  is  satisfied,  and  employing  large  argument  expansions,  one  may  show  that 


Ai(t)  -1/2 
■ Ai'(t)  ■ 


(1 


lt-V2 
4 ' 


). 


-IT  < arg  t < 0 


(25) 


Then  by  I.aplace  inversion  of  (24), 


_-ilT/4 


1 +• 


3/2  ilj 


60 


+ 0(y’/2) 


(26) 


Actually,  the  inversion  yields  for  the  factor  outside  the  braces  the  large  argument  approximation 
of  (i/2)  HQ(i)(ks)  the  Green's  function  for  an  infinite  perfectly  conducting  plane.  By  inserting  the  exact 
limiting  value  for  a — w,  the  formula  in  (26)  may  be  applied  as  well  in  the  near  field  of  the  source. 

To  trace  the  transition  from  concave  to  convex  curvature  when  the  distance  s along  the  surface 
remains  fixed,  we  allow  a to  change  continuously  from  positive  to  negative  values  via  the  complex 
excursion  Q < arg  a < ir  , To  keep  a = a|0  -^'|  positive,  it  is  necessary  simultaneously  to  have 
arg  |((>  -0'|  = - arg  a . Moreover,  one  must  continuously  deform  the  integration  path  in  the  complex 
t-plane  to  keep  the  integral  convergent  (expdyt)  oscillatory)  when  y is  allowed  to  become  complex 
according  to  the  rules  stated  above.  These  conderations  lead  to  a straight  line  path  along  which 
argt  = -Tr/3,  2n/3,  when  arg  a = w.  Changing  variables  t = 4 exp(i2iT/3),  one  finds 


iks 
_e 

2irka 


2/3  w (4) 


(27) 


where  y is  again  given  by  (18)  and  the  path  C„  proceeds  along  the  real  4 axis.  The  expression  in  (27) 
is  the  known  result  for  the  field  on  the  surface  of  a convex  perfectly  conducting  cylinder.  By  contour 
deformation  about  the  singularities  4p  at  wj  (4  ) s 0,  one  derives  the  creeping  wave  series,  and  by 
expansion  analogous  to  (25)  and  Laplace  inversion  the  limiting  transition  as  y 0 (HASSERJIAN,  G. 
and  ISHIMARU,  A.,  1972).  Thus,  one  may  track  the  field  continuously  as  the  curvature  changes  from 
concave  to  convex  between  fixed  source  and  observation  points  on  the  surface. 

G.  Generalisation  to  Variable  Curvature 

Invoking  the  principle  of  locality,  which  characterises  high  frequency  propagation  and  diffraction, 
one  may  generalise  the  preceding  results  to  accommodate  concave  surfaces  with  slowly  variable  radius 
of  curvature  a(s).  The  results  are  as  follows.  For  the  integral  representation  in  (21): 


iks 


4n[f(0)f(s)] 


r ikDt  Ai(t) 

T72  aTw  ***  ’ 


(28) 


C 


where  the  source  point  is  located  at  s = 0,  s is  the  distance  between  the  source  and  observation  points 
measured  along  the  surface,  and 


f(s)  = ( 


ka(s) 

2 


1/3 

) 


D(s)  . i 


f «lu 

or  f^(4) 


(29) 
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RN  in  (22)  is  treated  in  the  same  manner  (BABICH,  V.  M.  and  BOLDYREV,  V.  S.  , 1972).  As  D(s)  — 0, 
with  8 small,  one  may  derive  the  generalization  of  (26), 


G 

o 


(1) 


(ks)  1 + 


-itr/4 


TT  (kD)^/^  + 


O(kD) 


(30) 


The  ray-optical  series  in  (10),  with  (8),  becomes  (BABICH,  V.  M.  and  BOLDYREV,  V.  S. , 1972) 


ikD 

n 


e 


D 

n 


(31) 


where  is  the  divergence  coefficient  incorporating  the  effects  of  surface  curvature. 


j a^TIjTfpcos  e 
°1  ^ V 2fji2-a(lj+i2)co8ej 

(31a) 

°2  = ''°«®2<^+^2'^^3> 


•<f~~  -J  ^*1^  2*  3'*’®'! ^2^*1^*  2 ®1  ®2  " 2®'2'^°®  ®2 


- 2f  jaj  (f  j+lg)  Sj 


- 2/ ji  ja2  cos  02  + aja2l3  cos  Oj  cos  02 


(31b) 


with  a^  the  surface  radius  of  curvature  and  6^  the  ray  angle  with  respect  to  surface  normal  of  the  i-th 
reflection.  Finally,  the  asymptotically  approximated  (see  (3))  whispering  gallery  mode  series  in  (13), 


°oM 


iks 

2 [f(0)f(8)]‘/^ 


M 


z 


ikDt 

m 


t 

m 


(32) 


It  should  be  noted  that  (28)  cannot  be  used  when  a(s)  changes  continuously  from  positive  to  negative 
values  unless  the  source  is  located  near  the  inflection  point  a = oo  , where  (30)  and  its  convex  counterpart 
applies.  The  failure  of  (28)  for  arbitrary  source  location  on  the  concave  portion  arises  from  the  fact 
that  the  function  f(s)  -*'  «>  as  s nears  the  inflection  point.  On  the  other  hand,  D remains  finite.  The 
coefficient  nullifies  the  finite  integral. 

H.  Numerical  Results 

Some  numerical  calculations  in  the  literature  illustrate  the  utility  of  the  various  representations 
for  the  Green's  function  Gq.  Wasylkiwskyj  (WASYLKIWSKYJ,  W. , 1975)  has  evaluated  the  field  on  a 
concave  circular  perfectly  conducting  surface  by  summing  the  whispering  gallery  mode  series  and  the 
continuous  spectrum  contribution  in  (15)  (Figure  6).  Although  he  employed  this  method  of  computation, 
the  interpretation  of  the  result,  as  stated  by  him,  is  in  fact  achieved  via  the  ray-optical  formulation 
since  the  oscillations  in  the  curve  are  due  to  interference  between  the  direct  ray  and  reflected  rays. 
Wasylkiwskyj  observes  that  the  mixed  representation  in  (14),  while  useful  for  interpretation,  offers  no 
numerical  advantage.  It  appears  nevertheless  that  this  aspect  would  warrant  further  study.  When  ka  is 
changed  the  corresponding  curve  for  |Gq|  may  be  obtained  from  that  in  Fig.  6 by  a scaling  law 
(WASYLKIWSKYJ,  W.  , 1975),  thus  making  a typical  curve  universal  for  other  cases. 

When  ka  and  (or)  Y ia  relatively  small,  it  is  appropriate  to  explore  the  representation  in  (10). 

This  was  done  by  Buldyrev  and  Lanin  (BULDYREV,  V.  M.  and  LANIN,  A.  I. , 1975)  for  a concave  surface 
of  parabolic  shape,  for  which  the  asymptotic  results  could  be  checked  against  the  exact  series  solution. 
The  calculation  utilized  tabulated  values  for  Gq  in  (28)  or  Rj^  in  (22)  (when  generalized  as  in  (28)  to 
variable  curvature),  and  the  ray  series  in  (31),  with  N determined  from  (23)  (y  there  is  replaced  by 
kD).  One  observes  from  the  typical  curve  in  Fig.  7 that  the  asymptotic  representation  in  (28),  when 
modified  to  account  for  the  different  boundary  condition  on  Uq,  does  not  adequately  predict  the  field  for 
kD  so  large  that  a ray-optical  term  can  be  separated  out.  As  kD  increases  further,  each  formula  with 
a lower  N value  than  permitted  is  less  accurate  than  that  with  the  correct  N value.  This  demonstration 
shows  indeed  the  important  role  played  by  the  geometrical  optics  field. 


m.  CONCAVE  BOUNDARY  WITH  SURFACE  IMPEDANCE  CONDITION 

To  address  the  problem  of  ground  wave  propagation,  we  assume  that  the  ground  characteristics  can 
be  taken  into  account  by  a constant  surface  impedance  Z.  For  the  circular  cylinder,  this  implies  the 
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boundary  condition 


|^  = ikZ'G  at 

3 p 


p = a , ~ 


where  ^ is  the  impedance  of  free  space.  By  proceeding  as  in  Section  II,  one  may  show  that  the  following 
modified  expressions  apply: 

Instead  of  (1 ), 


1 - 

i(TTka)^  [h|^  *^*(ka)  - iZ'H^^^’cka)]  [jj^  (ka)  - iZ' (ka)] 


Instead  of  (4),  the  traveling  wave  expansion  involves 


H'  *‘*(ka)  - iZ'H*‘’(ka) 

H'  ’^*{ka)  - iZ'H*^’(ka)  ' 

V V 

with  corresponding  modification  of  (5).  The  saddle  points  in  (7)  are  unchanged,  and  the  ray-optical 
series  (8)  becomes 


2 ikD 

•~/A  r~T'  SS  (r  ) cos  w n 

_ _it/4  12  V / m"  n n e 

o Vek  (cosw  +Z')^  Jd~ 


where  is  the  boundary  reflection  coefficient  for  ray  species  n. 


cos  w - Z* 
n cos  w + Z* 


cosw„  = sin  , D^=  2a  (n+1)  sin-lfj^ 


The  whispering  gallery  mode  series  is 


ip  I0  -(|>'|  f 1 

e J (ka)  J'  (ka)  - iZ'J  (ka)  ‘ 

p ' 3p  I p ' p ' iv 


where  the  summation  extends  over  the  xeros  v defined  by 

m ^ 


J»  (ka)  - iZ<J  (ka)  = 0 
m m 


The  remainder  term  (11)  becomes 


ika|d)-(h'|  sinw  N+1 . , , .N+1  , 

^ _ 1 r e ly  y I (r  )(-l)  coswdw 

IST  d ^H^‘^>(ka)  - iZ'H|,^’(ka)J  ^J^(ka)  - iZ'J^(ka)j 


where  r is  given  by  (35a),  The  analogue  of  (21)  becomes 
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2/3  ^iYt 

2iir^(ka)  ' ^ [w^(t)  + iaWgW]  [Ai'(t)  + ioAi(t)] 

where  at  = Z'  Under  the  same  conditions,  one  has  for  Rj^  in  (22), 

^ 2/3  (.i)N+2 ^N+1  ^iyt 

2iir^(ka)  ^ w^(t)  + iaWgCt)  Ai'(t)  + io:Ai(t) 

C 


(42) 


(43) 


Employing  the  Wronskian,  the  field  may  be  put  in  the  form  of  (24) 


G 

o 


iks 
_e 

2irka 


<T)'  J 


2/3 


Aijt) j 

Ai'(t)  + iaAi(t) 


(44) 


The  perturbation  expansion 


Ai(t) 

Ai'(t)  + iQrAi(t) 


ia 

rt 


“ ,1  . . 3 


(45) 


leads  from  (44)  to  the  representation 


G 

o 


- [l  - y‘/^  a -2ia^y  + (^  + ia^  ) /T 


+ 0(Y^ 


(46) 


thereby  generalizing  (26). 
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FIG.  3 Whispering  gallery  mode:  modal  field,  modal  rays  and  caustic  {v^  = ka  sinw^) 


FIG.  4 


Integration  path  and  singularities 
X --  zeros  Vp  of  H^^^^ka); 


in  complex  v -plane 
. --  zeros  Vjjj  of  JJ,  (ka) 


5 Integration  paths  and  singularities  in  complex  w plane  {v  = ka  sinw).  C corresponds  to  C in 
Fig.  4 and  can  be  changed  to  C because  of  synunetries  in  the  integrand  in  (5).  C'  in  turn  can 
be  deformed  into  the  steepest  descent  path  through  w,„.  The  pole  singularities  Wm  arising 
from  the  zeros  of  (ka)  do  not  affect  the  asymptotic  evaluation  of  (5).  The  poles  Wp  are  not 
shown. 
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DISCUSSION 


K. l).  Becker : In  1051,  B.  Priedman  showed  that,  in  the  case  of  an  ideally  conducting  sphere, 
the  continuous  spectrum  is  vanishing.  Is  there  in  your  case  a connection  with  the  sign 

of  curvature  of  ai  ideally  conducting  surface  and  the  existence  of  a continuous  spectrum? 

L,  B, Felsen : The  existence  or  not  of  a continuous  spectrum  depends  on  the  choice  of  the 
representative  theorem.  In  the  problem  which  I discussed  in  cylindrical  geometry,  the 
field  was  represented  in  terms  of  angularly  propagating  waves  in  an  infinity-extended 
angular  space.  Therefore,  the  spectral 'theorem  has  to  he  developed  for  the  radial  domain. 
For  the  interior  of  the  cylinder,  with  the  origin  (cylinder  axis)  accessible,  one  finds 
both,  a discrete  and  a continuous  spectrum.  When  the  origin  is  shielded  by  another  cy- 
lindrical boundary,  the  spectrum  is  purely  discrete. 


